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Abstract. Doppler refiectometry spatial and wavenumber resolution is analyzed 
within the framework of the linear Born approximation in slab plasma model. Explicit 
expression for its signal backscattering spectrum is obtained in terms of wavenumber 
and frequency spectra of turbulence which is assumed to be radially statistically 
inhomogeneous. Scattering efficiency for both back and forward scattering (in radial 
direction) is introduced and shown to be inverse proportional to the square of radial 
wavenumber of the probing wave at the fluctuation location thus making the spatial 
resolution of diagnostics sensitive to density profile. It is shown that in case of forward 
scattering additional localization can be provided by the antenna diagram. It is 
demonstrated that in case of backscattering the spatial resolution can be better if 
the turbulence spectrum at high radial wavenumbers is suppressed. The improvement 
of Doppler refiectometry data localization by probing beam focusing onto the cut-off is 
proposed and described. The possibility of Doppler refiectometry data interpretation 
based on the obtained expressions is shown. 
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1. Introduction 

Plasma rotation velocity measurements are of great importance for understanding 
transition to improved confinement in tokamaks and physics of transport barriers. 
Extensively used nowadays for such investigations is Doppler reflectometry EJ [Hj. 
This technique provides measuring fluctuations propagation poloidal velocity which is 
often shown to be dominated hj E x B velocity of plasma [3]. In using this method a 
probing microwave beam is launched into the plasma with finite tilt angle with respect 
to density gradient. A back-scattered signal with frequency differing from the probing 
one is registered by a nearby standing or the same antenna. The information on 
plasma poloidal rotation is obtained in this technique from the frequency shift of the 
backscattering spectrum which is supposed to originate from the Doppler effect due 
to the fluctuation rotation. Spatial distribution of the scattering phenomena, which 
is usually assumed to occur in the cut-off vicinity, is the key issue for the diagnostic 
applications. 

Some numerical simulations undertaken (see j3] and references there, ^) and 
analytical results jH] demonstrate the possibility of the measurements localization by 
the cut-off but the problem of Doppler reflectometry locality remains still open and 
lacks comprehensive analytical treatment. This stems from the fact that to be sure that 
main contribution to the registered signal is made by the cut-off vicinity one should 
compare it with integral contributions of distant from the cut-off regions which can be 
substantial depending on density profile. 

This paper attempts to clarify this problem. Here we present a theoretical 
investigation of the Doppler reflectometry signal dependence on the turbulence 
distribution with respect to the cut-off. We consider slab two-dimensional model 
(see figure C]). This simplification allows us to perform straightforward analytical 
treatment and obtain explicit reliable expressions for the scattered signal which can be 
easily used for estimation of diagnostics locality and experimental data interpretation 
without complicated and time consuming numerical calculation using full-wave codes, 
etc. The model considered can be readily applied to large in comparison with 
probing beam width, vertically elongated (ITER-like) plasma. The limitations of this 
approximation are discussed below. 

The paper is organized as follows. In section |21 the consideration is carried on 
in the geometrical optics approximation for arbitrary plasma density profile. This 
approximation fails to hold in the cut-off region. So in this region the analysis is made 
assuming the density profile to be linear and using exact expressions for the probing 
and scattered wave electric field given by Airy functions (sectional). Some numerical 
examples are given in section |3] Brief discussion of the model considered and results 
obtained is provided in section El Finally, a conclusion follows in section IHl 
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Figure 1. Diagnostics scheme. 1 — emitting and receiving antenna, 2 — additional 
receiving antenna (to be discussed in section I^J. 



2. Reflectometry signal in WKB— approximation 

In this section the reflectometry signal is obtained in geometrical optics (or WKB) 
approximation. The plasma is assumed to be nonuniform in x (radial) direction and 
uniform in y (poloidal) direction. External magnetic fleld is supposed to be along z axis. 
0-mode Doppler reflectometry is considered but the flnal results can be easily adapted 
to X-mode reflectometry applying corresponding expression for wave radial wavenumber 
dependence on the radial coordinate. 

A received signal is obtained using reciprocity theorem and is assumed to be 
created by single scattering (linear) mechanism. The frameworks of this approximation 
and experimental means to check its applicability are discussed in section El The 
scattering signal averaging is made taking into account radial statistical inhomogeneity 
of the turbulence. Then a scattering efficiency is introduced and analyzed. 

2.1. Scattering signal 

We consider normalized antenna electric fleld in the following form 



where factor f{ky) takes into account the antenna pattern describing antenna radiation 
in vacuum 



The vacuum antenna fleld Eq differs from Ea by the absence of the reflected wave 
contribution. 

Here and further we suppose the probing wave to oscillate at frequency uj and omit 
the corresponding term in equations. Radial electric field distribution is described by 
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the function W{x, ky) which is determined by 

W" + kl{X,ky)W = Q (1) 

where the square of ordinary wave wavevector radial projection is given by k1{x, ky) = 
k'^{x) — ky = [cu^ — ujp^{x)]/c^ — ky. Thus ordinary wave electric field in WKB- 
approximation has the following form [7j: 

W{x, ky) 



C2 ^ 
'^2:^ 7 y . 



X cos 



exp 



Xc{ky) 



ITT 



kx{x , ky) dx — 



where the turning point Xc{ky) is determined by the equation 

kx [^c{ky), ky] = 

and X = corresponds to the plasma border. 

Using reciprocity theorem [HIIZI, an amplitude of the received scattering signal at 
frequency Ug = u + Q, where Q denotes the frequency of the fluctuation caused the 
scattering, can be written as 

o2 



Asi^s) 



le 



6nQ{f)E^{r) df 



Here Pi is the probing wave power. Following the model considered the fluctuations are 
assumed to be long enough along the magnetic field direction z, so their dependence on 
z can be neglected. Introducing the density perturbation Fourier harmonic 

"+00 



(5n(x, q, Q) 



6n{x, y, fi)e-'"^-^'?Mxd?/ 



we obtain 



As{uJs) 



me 
2mpU 



+00 



dky dxdq 
(27r)3 



5n{K, q, n)f{ky)f{-ky - q)C{x, q, ky 



(2) 



The efficiency of scattering off the fluctuation with radial wavenumber x and poloidal 
wavenumber q has the following form [7j 



C(x, q, ky) 



W{x, ky)W{x, -ky - g)e^"^da; 
dx 



{m+n)n / A 



^/kxix, ky)k,jr.{x, ky q) ^^^^^^ 
= XX + m(f){x, ky) + n(f){x, —ky — q) + 0(0, ky) + 0(0, —ky — q) 



{x, ky) 



Xc{ky) 



k^i^x , ky) dx 



The scattering efficiency C{>i,q,ky) is an integral of the oscillating function. Main 
contribution to this integral is made by stationary phase points. Following j7] to estimate 
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this integral let us calculate this contribution. Equalizing the derivative of the phase to 
zero one obtains the Bragg condition 



X ■ 



In nonuniform plasma this determines the scattering point 



X q{2ky + q) 
2 2x 



1 2 



The figures m, n are related to different cases of scattering. 



m 



sgn 



x2 + q{2ky + q) 



n 



X 



sgn 



x^ - q{2ky + q) 



X 



(3) 



(4) 



(5) 



As it will be shown below case m = n corresponds to the backscattering (BS) and 
m = —n accords to the forward scattering (FS). The sign of m is related to the 
scattering before (m = —1) and after (m = 1) the cut-off in respect of the probing 
wave propagation. Equations (0) express the fact that one fluctuation can scatter the 
wave only once. The situation, for example, when the wave with fixed ky is forward 
scattered off the fluctuation far from the cut-off and the same wave can be scattered 
backward off the same fluctuation near the cut-off is impossible. This circumstance 
will allow us below to separate the contributions of forward and backward scattering 
processes. 

The final expression for the scattering efficiency in WKB approximation takes the 

form 



where 



C(x, g, ky) = 2W-— ^exp 



[dk'^{x)/dx\ 




i^rnn{,x^) + — sgn(mnx) 



is local Airy scale. Introducing local density 
variation scale = [dhine^x) /dx\x=x,Y^ we obtain for ordinary probing wave = 

On substituting the obtained expression for the scattering efficiency into I 
obtain the received signal 



we 



+ 00 



dky dx dg 
(27r)3 



.3/2 



6n{x,q,Q)f{ky)f{-ky - q) 



X ■ 



10 



(6) 



where o > determines how the integration path goes around the singularity. Here 



\E' (x^ , ky) 



\E'mn, where m,n is determined by ©• 



2.2. Scattering signal analysis: integration over ky 

We calculate the integral over ky in (jH)) using saddle point method. It is especially 
efficient if the cut-off is far enough from the antenna and does not coincide with the 
focal point of the antenna beam, so that the ray tracing consideration is applicable. 
In this case the saddle point k* is determined by the stationary point of the phase 
$ = ^(x*, ky) + Bigf{ky) + axgf{—ky — q) and the \f{ky)f{—ky — q)\ dependence on 
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ky is negligible. Corresponding criterion to distinguish this and opposite case will be 
formulated below. 

In this "ray tracing" case the stationary phase condition 

d$ [x^{ky), ky] 



dky 



(7) 



can be easily interpreted. At first we consider the meaning of /dky term. Recollecting 
that 



and 



one obtains 



k^{x^ ky^ — k (x) ky (8) 



d(f){x, ky) _ dx' 



''^^y / n — TT ~ '^y^y + 



dky J kx{x^ky) kx{x, ky q) 

Jo kx{x, ky) Jq kx{x, —ky — q) 

Taking into account that (jH)) yields the following relation between the projections of 
group velocities Vg = du/dk of incident and scattered waves: 



h 7/"*) ~k 

•Jgy^ _ rvy Ugy _ ft^ 

kx{x,ky) -y^*) kx{x,—ky — q) 



•Jgx --^y-i -y/ ugx 

we obtain 

rxciky) (i) rxc(ky) ii) rxc{-ky-q) (s) ^xd-ky-q) (s) 

Jx, Vifi Jo Vli Jx, V'gj Jo 

Now one can see that Ay = —d'^ /dky corresponds to shift of the ray trajectory along 
y direction when it returns to the antenna. In less general case it was mentioned in [S]. 
For instance, if m = —n = 1 then (fTIH) takes the form 

Jo Vgx Jxc(ky) -Vgx Jx, VgJ 

and {—d'^/dky) is the shift of the ray trajectory in case of forward scattering after the 
reflection off the turning point. Here the meaning of the figures m,n ^ announced 
above becomes clear. 

Accounting for the influence of the wavefront curvature at the antenna given by 
axgf{ky) 7^ we consider gaussian antenna pattern 



f{ky) = ^20Fpe-('' /2 (^^^ 

In case of 

cTZ/uo > p2 (12) 
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parameter TZ has a meaning of a wavefront curvature radius at the antenna. In (lllj) /C 
corresponds to the antenna tilt (/C = cj/csin??, where denotes tilt angle in respect of 
the density gradient). This allows us to obtain condition in the form 
d$ cTZ 

_._A, + -(2l, + ,) = (13) 

Taking into account that the ray with poloidal wavevector component ky is radiated 
from the position y{ky) = —kyTZc/ut of the curved wavefront at the antenna which is 
assumed to be situated in the axes origin, it is easy to show that (fT^ determines a 
ray trajectory which returns to the proper point of the antenna wavefront after the 
scattering off the fluctuation with wavevector (x, g). 

To get explicit expression for the stationary point position k* we use paraxial 
approximation. Supposition of small divergence of the antenna beam \ky — IC\ ^ cu/c 
which holds usually in the experiments on Doppler reflectometry allows us to write 

k^-IC^ 
kx{x, ky) ^ kx{x, K,) — 2^ 

and neglect the dependence on ky in denominators in (jH)). This allows us to obtain 

1^ = [(^ + n)A{x,) + 2Ao] - — [nA{x,) + Aq] 

OKy UJ U! 

r-Xc{/C) 



Using (|7j) we get the stationary phase point position 



, -q/2, m=n 

K = { (14) 

^ ' -g/2[l+r^A(x,)/(Ao-7^)], m=-n 

for the BS and FS respectively. 

In the opposite case when ray trajectories consideration is not valid k* is determined 
by the antenna pattern amplitude 

d\f{ky)f{-ky-q)\ 



dky 







yielding k* = —q/2 for an arbitrary antenna pattern. 

The assumption of gaussian antenna pattern (fTT|l allows us to evaluate the saddle 
point position in general case 

_ q Aq - 7^ + nA{x^) - iP 
y~ 2Ao-7^+ (m + n)/2-A(x,) -iP 

where V = ujp^ /c and formulate the criterion in question. The ray tracing consideration 
is valid if <^ c|Ao — TZ\/uj and thus the focal point is not too close to the cut-off. If 
the focal point is situated in the cut-off ^ c|Ao — lZ\/uj we come to the expression 
K = -q/2. 
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Using stationary k* calculated above and performing the integration we obtain the 
scattering signal in the following form 



A, = 27r- 



rUpC^ 



+00 



dxdq 
"(2^ 



5n(x, q, fi) 



/ q)] f [-k;{x, q) - q] tj^^^^^ 



(x.,A:*) 



where 



A(x,g) 



~ ^ 2 dkl 



-1/2 



2.3. Scattering signal averaging 

We consider the turbulence to be slightly inhomogeneous along x direction so that the 
density fluctuation correlation function takes the form 



{5n{x)5n{x')) = br? 



X + x' 



dx 
27r 



n [ x,q, fi. 



X + x' 



^i>t{x—x') 



(15) 



This representation is applicable when the turbulence correlation length along x axis 
icx is much smaller than the turbulence inhomogeneity scale. It allows us to take into 
account the dependence of turbulence on the radial coordinate and still describe it 
using wavenumber spectrum \n[>i,q,Q,{x + x')/2]f. Supposing the turbulence to be 
stationary and homogeneous in ?/-direction and using we can represent correlation 
function of spectral density in the form 



((5n(x, q, f2)(5n*(x', q', Q')) 



dx6n (x) 



n 



X + x' 



g, X 



where the integration is performed over all plasma volume. It allows spectral power 
density of the received signal to be represented in the following form 



p{ujs) = {AsAs) = Pi 



dx 5n'^{x)S{x) 



(16) 



where Ag is complex conjugate to Ag. 

The scattering efficiency S{x) introduced here takes the form 



S{x) = 471^ 



dxdx' dq 
(27r)3 



n 



X + V 



g, VL, X 



X 



/ [fc;(x, q)] f [-k;{x, q)-q]f [k;{>,\ q)] f [-fc;(x', q) - q] 
V xx' 

and determines actually the Doppler reflectometry locality and wavenumber resolution. 
Performing the integration over x — V using stationary phase method we obtain 

.2 N 2 



S{x) ^ (27r)3/2 



dxdq |_ 



n{x,q,n,x)\ |/ [A;*(x,g)] I |/ [-A;*(x, g) - g] 
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X ii ; ; exp { 'J^^ - -sgnx:^(x, q) \ (17) 



klv^^J I 2x;^(^,g) 4 
where x'^^(x, g) = dx^{x,q)/dx and the scattering point position x*(x, g) is given 
by ©. 

One can see that due to the oscillating term exp {i [x — a;*(x, q)f / (2x'^^(x, g))} the 
main input to the integral (fTTj) is provided by the stationary point x = x^. It corresponds 
to the contribution of the fluctuations producing scattering just in the point x. 

Here we suppose the fluctuation spectral density |n(x, g, f2, x)^ to vary with k slow 
enough in comparison with the oscillating term. In general case four terms arise from 
condition x^{x, q) = x which is equivalent to 

k[x,{x\q)\ = k{x) (18) 

where k{x^) is determined by (jlj) and ky = k* is substituted according to section 1221 
Solving (fTH|) and taking interest in "stationary" x* we obtain 



>^*m,n = mk^ix, k*y) + kl{x, k^) - q{2k*y + q) (19) 

Figures m and n have the same meaning as above but here and further they are 
independent parameters and specify corresponding solutions of (fTH)) . 

To distinguish contributions of different scattering types we represent the scattering 
efficiency in the following form 

^(^) = ^ (;;^) ^ [Sbs{x, q) + Sfs{x, g)] (20) 

and discuss the properties of back and forward scattering efficiency separately. 

2.4- BS efficiency 

To obtain the two terms corresponding to the backscattering which occur before or after 
reflection of the probing wave we substitute m = n and k* = — g/2 to (fT^ according 
to jUl). This yields 



2mA;.,; ( x, -| ) , m = ±1 



and 



^*kx (^x, 2^ 



.2 ^ -1/2 



I q)\' = \p' + ^[Ao-n + mAix)f 

It allows us to obtain BS contribution in question 

Q l/(~g/2)l^ \^ \h[2mkx {x,K:) ,q,n,x]f 

i^x Kx, K.) ^^^^ + c2 [Ao - 7^ + mA{x)f /uo^ 

It can be seen that the contributions of the BS before (m=n=-l) and after (m=n=l) 
the cut-off have different amplitudes due to the diffraction effect dependence on length 
of the ray trajectory from emitting to receiving antenna. 
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We consider main features of the BS efficiency (|21|) obtained. The first factor 
determining the locahzation of BS signal which appears to be typical for the fiuctuation 
refiectometry jHJ is k~'^{x). Just recently [9J the same factor was revealed in analysis of 
correlation matrix of phase perturbations of fiuctuation refiectometry signal in nonlinear 
regime. It is maximal in the vicinity of the cut-off where in the WKB approximation 
it has singularity. This singularity saturation will be treated in details in section El 
This maximum explained by the growth of the probing and scattered wave electric filed 
in the cut-off provides the technique with spatial localization, however, the decay of 
when leaving cut-off is not fast enough to guarantee suppression of the signal 
coming from wide edge region. Moreover, for quite a few density profiles (e.g. linear 
and bent down ones) the integral of k~'^{x) over x does not converge due to far from the 
cut-off regions. So plasma periphery contribution can be essential in these cases. The 
illustration of this effect is given in section 0] 

Additional possibility to increase selectively the BS signal coming from the cut- 
off can be provided by focusing the beam of probing and receiving antennae to the 
cut-off region. This corresponds to 7^ = Aq in (j^H). In this case the focusing causes 
additional growth at the cut-off of probing and scattered wave amplitudes leading to the 
scattered signal enhancement similar to that predicted for the cut-off or upper hybrid 
resonance [TU] . 

The two mentioned effects increasing the probing wave electric field substantially in 
the cut-off should enhance the locality of the Doppler reflectometry diagnostics making 
it less sensitive to the backscattering in the edge plasma. 

The backscattering locality can be also better for some fluctuation radial 
wavenumber spectra. One can see from PHI that if the short-scale fluctuations are 
suppressed enough in the spectrum the backscattering contribution will be essential in 
the cut-off vicinity only. 



2.5. FS efficiency 



To calculate the FS terms we use ray-tracing consideration. Substituting m = —n 
to |T9|l we obtain 

g^A(x) 



' 2K{x,k;){Ao-n) 

It is noteworthy that both xj^ _^ have the same sign and differ only slightly due to 
different values of k* only. 
Finally we obtain 



x^ 



2c(Ao - 7^) 



2 



k^i^x, ky)kx{x, ky q) 

c\A,-nf^-"' 



ujk{x) 



2cL^:kx{x, k*y) 



1 



2cLif:kx{x, ky) 



uA{x) 
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Thus the sum of two similar forward scattering terms corresponding to the 
scattering before and after the turning point takes a form 



Spsix^q) 



p4 + c2(Ao - ny/io^ 
A{x) 



f 



An -7^ 



/ 



1 



An -7^ 



1 



n 



2k{x){Ao-n) 



(22) 



f^x (^) ^) 

The main reason for the mentioned similarity is equal trajectory length in this case. It 
is worth noting that if 7^ < Aq the forward scattering efficiency Sps is contributed to 
by the fluctuations with positive radial wavenumbers only. 

Substituting gaussian antenna power diagram (jllj) into (j22j) we obtain expression 
for the factor describing antenna pattern influence in the following form 



/ 




A{x) 



iq + 2icy + 



qA{x) 


2" 




_Ao-7^. 




} 



(23) 



Additional contribution to the forward scattering efficiency can be provided by 
singularities of the term |A(x, g)|^ / |x| in ()17|1 . As it is shown in [Appendix A| it is 
small at substantial distance from the cut-off where the condition 



ck{x) 



UJ 



> P 



12 



a 



1/2 



c(Ao - 7^) 
\Ll/ne{xc) ■ d'^ne{x)/dx'^ 



(24) 

characterizes the nonlinearity of the 



holds. Here a 

density profile. If the cut-off is not close to the focal point, so that p^yuj/[c{Ao — TZ)] -C 1 
the expression (j^^ for the FS contribution appears to be valid where nonlinear 
corrections to k'^{x) are still small. Closer to the cut-off where inequality is not 
satisfied the contribution of the branching point to the integral in (fT7j) is important. 
It is calculated below taking into account that the density profile in this region can be 
treated as linear. The corresponding contribution doubles the result for the FS signal. 

Similar to ()21|) FS efficiency ()22j) is also proportional to k~'^{x). However unlike ()2ip 
the cut-off contribution is not enhanced there by probing wave focusing. In this case 
the focusing merely compensates the refraction of probing beam leading to FS signal 
amplitude growth in all plasma volume. It can be seen from (j^^ that FS signal is 
contributed to by long-scale fluctuations, which disable the third possible localizing 
factor discussed for BS — turbulence spectrum. 

Nevertheless, according to ((221), the extra localization of the forward scattering 
can be due to the fact that FS signal coming from plasma volume is received in Doppler 
reflectometry is received by the antenna pattern periphery. Supposing gaussian antenna 
beam to be wide enough, the integration over q yields the following estimation of the 
forward scattering efficiency 

„2 \ 2 p+oo 



Sfs{x) 



dq 
2^ 



SFs{x,q) 
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X 



n 



' 2 {271 f^ujp 
c(Ao-7^)A;2(x 

2}C^A{x) 



exp 



pJCAjx) 



— 2/C, r2, X 



[kix){Ao-ny 

Thus the factor describing the antenna pattern effect decreases rapidly in moving off 
the cut-off under condition the beam is wide enough or sufficiently tilted. If the antenna 
beam focusing to the cut-off is provided this localizing factor can be even more essential. 



3. Scattering signal in the cut-off vicinity 



The explicit expressions for the scattering efficiency (j2^, were obtained in the 
previous section using WKB and ray tracing consideration. These expressions possesses 
singularities when validity conditions for WKB approach 
dfc(x) 



k-\x)- 



dx 



< 1 



(25) 



and ray tracing consideration 
c(Ao - 7^) 



UJ 



are violated in the cut-off vicinity. 

To analyze the scattering efficiency in this region a more rigorous approach is 
needed. To do that we recollect that in the cut-off vicinity the density profile can 

1 /3 

be represented as linear k'^{x) = {xc — x)/i^ where £ = {c^L/u'^) is the Airy scale and 
L = [d\nne{x) / dx\x=xc]~^ is local density variation scale in the cut-off position. In this 
case the criterion ()25p takes the form k{x)i ^ 1. 

Then in the cut-off vicinity the radial distribution of the ordinary wave electric field 
has the following form [7j: 



W{x, ky 



— exp 



c 



Xc{ky) 



k^i^x 5 ky^ dx 



ITC 



X 



exp 



ip 



+ + k'yHP 



dp 



where ^ = {x — Xc)/i- 

According to the reciprocity theorem jH] the scattering signal has the following 
structure 

2 ^+°° dfcj^dxdg 



A 



le 

ArripU! 



dx 



{2ny 



■ 6n{K, q, ^l)e 



\>c{x-Xc) 



X f{ky)f{-ky - q)W{X, ky)W{X, -ky - q) (26) 

As we have already done assuming moderate spatial inhomogeneity of the 
turbulence (fT^ and calculating the integral over ky in (j2Sl) by stationary phase method 
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we transform the registered signal spectral power density to the form For the 

scattering efficiency we get 



S{x) 



y + c2(Ao -7^)2/cJ2 yrrieC 

\n[>£, q, iZ, x)\ 



R{x, X, q) 



oo 
oo 



(271) 



/(-f 



R{x, X, q) 



dr 



2 / , . x2 



exp 



(27) 
(28) 



is: = 2lkr X 



where the following notation is used 

/9- p ft - ^^""j" ^° - ^ - 

^ ' u; (Ao-7^)2 + p2' ^ Ao-7^' 

The expression for R{x, x, g) can be simplified in two cases. In the first one K ^ 6 
and the position of the turbulence x is far from the cut-off. Omitting oscillating terms 
hardly having an effect in integrating over f3 we can distinguish two characteristic ranges 
of fluctuation radial wavenumbers. The fist one corresponds to <^ K. Considering 
such fluctuations, which are responsible for the forward scattering we can approximate 

'1 



R ^ R^^\x,x,q) = 27rJc 



-K' \/3-e\ 



eX2 



(29) 



The second group of wavenumbers corresponds to |/?| ~ These fluctuations provide 
the backscattering in this region. Corresponding expression of the scattering efficiency 
has the following form 



24/3„ 

R^R^^\x,x,q) = —r-Ai 



22/3 



(30) 



Assuming the spectral density |fi(x, q, Q, x) \ to vary with x slow enough and neglecting 
this variation when performing the integration over x which is correct for distances 
from the cut-off Xc — x > l^x where lex is radial correlation length of the turbulence, 
we represent the scattering efficiency in form ()2Up and get the forward scattering 
contribution 



Spsix.q) = 4 



-1/2 



2{pqLf 



X exp 
and backscattering one 



l/(-g/2)r 

kl (x, -q/2) 
ck^ (x, -q/2) 



Lcq^ 



An-7^ 



uj (Ao-7^)2 + p2 



,q,VL,x 



SBs{x,q) 



(Ao - 7^)2 + p2 L ^ 

|/(-g/2)r 



(31) 



kUx,-q/2) 



Y,\n[±2k,{x,-q/2),q,n,x]f (32) 



This BS efficiency expression matches when leaving the cut-off corresponding WKB 
formula (PT|) obtained for the arbitrary density profile whereas the FS contribution ()3ip 
exceeds corresponding WKB result (|22|) by the factor of 2. Taking into account the 
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discussion of the previous section we can write approximate formula describing the 
transition from ()22|) to (jHT|l which happens when nonhnear correction to the density 
profile decreases. 

1 [pqA{x)f 



1 + 272 



X 



4 c\Ao-nf 



q_ 

'2 

-1/2 



n 



'q^K^x) Ao-7^ 



_2k{x) (Ao-7^)2 + p2 
ck{x) /{up) ^ac{Ao-n)/{12uj). 



q, fl, X 



(33) 



where as above V = up'^/c and 7 

Expressions (jSH), (jHS)) and for the scattering efficiency describe the transition 
from the ray tracing consideration to the case of probing beam focusing to the cut-off 
both for the gaussian and arbitrary antenna pattern. In the last case we determine 
parameters p, IZ as: 



P 



-3? 



r(/c) 



n 



'Pilcy c/2(/c) 

In case when ray tracing approximation fails to hold, which can be provided by large 
antenna beam width > cAq/u or by the focusing to the cut-off 71 ~ Aq, the FS 
contribution takes the following form 

2" 



Spsix^q) ^ 4 [pk{x)] 



-2 



/(-I 



exp 



1 fcq\'' 



2 \uj/ 



A(x) 
P 



\n [0, q, fl, x] 



Equation (|7T|) for BS contribution holds true in general situation (V ^ \Aq — TZ\) 
for arbitrary antenna pattern providing the redefining of p and TZ mentioned above is 
made. 

The second important case to be considered is when the turbulence is situated near 
the cut-off K < 1 where WKB scattering efficiency has singularity. Main contribution 
to the scattering in this location can be shown to be due to fluctuations with \f3\ < K. 
It is given by 



X, q) ~ 27rJo 



2 

°°dr 
r 



(ir' - \i3 - e\ 



2 — 



sm 



-e 



-/32)/2 



TT 



This expression reveals a scattering efficiency maximum to be situated near K ^ 2.3 
which corresponds to k{x) ~ 1.15/£ or Xc~-x ~ 1.3i. A numerical calculation of function 
R{x, X, q) with following integration over x confirms this result. Figure |21 represents 
numerical results (black dots) in case of short-scale fluctuations icx ^ ^- Here for the 
sake of clarity the fluctuations spectral density is supposed to be constant in ^I7\i during 
the integration over x. 

Additionally this computation shows WKB formulae (PT|) . for scattering 

efficiency to be valid up to the maximum providing k^ in denominators is replaced 
by 



k^ - 0.5/£ 



(34) 
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Figure 2. Scattering efRciency calculated numerically (•) and analytically (solid line) 
vs. K = 2k^{x,-q/2)L 



(solid line in figure Ej). Oscillations near this solid line correspond to the oscillating 
term omitted in (j^^ which can be shown not to contribute essentially to the resulting 
scattered signal magnitude due to averaging provided by slow spatial variation of the 
fluctuation amplitude. Fast decay of S for kx < 0.5/ i is caused by the probing wave 
field decrease in the evanescent region. 

In integration over x in expressions (fTTjl . (f^ . (jHn|l we neglected the influence of 
the fluctuation spectral density \n{>i, q, Q, x)f. Now we consider the case £cx ^ ^ when 
this approximation is not valid. In this case of long-scale turbulence the backscattering 
contribution is small in comparison with forward scattering one and can be neglected 



S{x) 



/ 7r^Ps{x,q) 

J-oo ^TT 



According to ^T7\i . the scattering efficiency takes the form 

2 



SFs{x,q) 



2c 



f - 



dp 



n 



q, Q, X 



K^\/3-9\ 



(35) 



The characteristic scale of the spectral density variation with (3 = >ci is i/icx- The 
corresponding scale for the term Jq[1/2K'^ \P — 9\] is 6P ~ K'"^. In evaluating last 
integral in two cases can be distinguished. Far from the cut-off at Xc — x ^ £cx the 
integral converges due to the Bessel function at /5 ~ 5/3 -C £/£ca;- It gives fc~^-behavior 
for the scattering efficiency which was obtained above ()31|) . 

In the opposite case Xc — x< icx the integral in (jHKjl converges due to the turbulence 
spectrum at /? ^ i/icx leading to the saturation of the singularity k~'^. This saturation 
is described analytically in the case of gaussian fluctuation spectral density 



\n (x, q, Q)f = 2y/TTicx \n {q, Q) f e 
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Figure 3. Approximation of scattering efRciency for £cx ^ i (solid line £cx — 2 cm, 
£ = 0.5 cm) and £cx < ^ (dashed line) vs. K = 2kx{x, ~q/2)£. 



which allows the integral over (3 in fl35p to be calculated exactly 

2 



d/3 



n 



2 / ^ 

27r£ |n (g, f2, x)| exp 



4,^2 \ ( K^oi 
lo 



where Iq is modified Bessel function. 

The difference in scattering efficiency behavior for small-scale [icx <C ^) and long- 
scale turbulence [icx ^ ^) is illustrated by figure El In our case the transition to the 
plateau occurs in the point K'^ ~ ^icxl^ which gives 

Xc X ~ ^cx 

To summarize, in the case lex -C ^ in approaching the cut-off the scattering efficiency 
grows, has maximum in the point Xc — x ~ 1.3£ and diminishes when Xc — x < £ due 
to probing wave field decrease in the evanescent region. When i^x ^ I the scattering 
efficiency has a plateau in the region Xc — x < i^x- 



4. Estimation of experiment locality 

On performing the correction ()34j) the WKB- formulae ()21|). ()22|) can be used up 
to the cut-off vicinity. We consider simple model illustrating main properties of 
the scattering efficiency obtained. Geometrical parameters taken correspond to Tore 
Supra experiments T|, where cut-off being situated in the antenna near- field zone 
(cj/c ~ 12.6 cm~^, p ~ 14 cm, distance to the cut-off L ~ 20 cm). Let us assume bent 
down plasma density profile (figure EJa)) similar to observed in Tore Supra |n| and 
step-like plasma poloidal velocity distribution (figure IH^b)). For the sake of simphcity 
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Figure 4. Signal spectrum evolution, (a) Assumed density profile. Horizontal axe is 
equivalent to the major radius. Probing is performed from right-hand side, (b) Poloidal 
velocity profile (solid line), cut-off positions (•) and velocity estimated using frequency 
spectrum shift (o). (c) Signal spectrum corresponding to the cut-off positions marked 
on (b). 

we assume the turbulence level to be uniform (which looks like real situation when 
relative turbulence amplitude 6n/n increases in approaching plasma periphery) and the 
wavenumber spectra to be gaussian. 

The probing is performed at different frequencies and therefore with different cut- 
off positions. We assume probing antenna to provide constant tilt angle i) = 11.5° and 
will calculate corresponding probing wave poloidal wavenumber K, = uj/c sin d for each 
frequency. 

We take into account scattered signal frequency shift due to the Doppler effect 
in the final expression for the scattering efficiency therefore assuming arising additional 
dependence on radial coordinate \h[. ...... ,VL — qv{x)] |^ to be slow enough. We perform 

the integration over q in (jSH), (jHH)) assuming the antenna pattern (fTT|) to be wide enough 
to determine the behavior of integrand and calculate spectral power density (fT6|l of the 
registered signal. Spectra obtained are represented in figure Efc). 

Turbulence correlation length was taken small Cc ~ 0.1 cm, Airy length being 
i ~ 0.5 cm. So the turbulence spectrum did not improve enough the locality and 
one can see that contribution of long area with v = —0.5 ■ 10^ cm/s (see figure |3I^b)) 
is essential for spectra 1-3 (figure IH^c)). It can be seen that these spectra Doppler 
shifts do not accord the velocity in corresponding cut-off positions (see figure Efd)). 
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Figure 5. Signal spectrum evolution, (a) Assumed density profile. Horizontal axe is 
equivalent to the major radius. Probing is performed from right-hand side, (b) Poloidal 
velocity profile (solid line), cut-off positions (•) and velocity estimated using frequency 
spectrum shift (o). (c) Signal spectrum corresponding to the cut-off positions marked 
on (b). 

It is necessary to move cut-off deep inside the plasma (spectra 4, 5) to provide the 
dominance of the region behind the velocity step point. Thus poor locality of Doppler 
reflectometry associated with bend down density profile in this situation can obscure 
the velocity distribution. 

Diagnostics localization dependence on plasma density profile is also demonstrated 
by the next example simulating a transport barrier. We consider plasma density 
profile (figure Efa)) bent down in plasma periphery and bent up in the core. The 
poloidal velocity profile has high gradient in the "barrier" region (figure Efb)). One can 
see that up to cut-off position 3 frequency spectrum shift (see spectra 1-3, figure Uljc)) 
accords to negative poloidal rotation velocity corresponding to the plasma periphery. 
But on crossing the twist point by the cut-off position the signal grows and the frequency 
shift changes to that corresponding to the local value of the velocity in the cut-off region 
(cf. spectra 4, 5, figure El^c) and figure Efd)). This illustrates the better locality of the 
method when used on bent up density profile. 

To illustrate antenna focusing influence we consider the density profile of DIII-D 
tokamak plasma with internal transport barrier (figure IHf a)) |T2]. Here we take into 
account the distance between antenna and the plasma, which was assumed to be 
equal 1 m, and suppose that the probing is performed with narrow antenna beam 
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Figure 6. Antenna focusing influence, (a) DIII-D density proflle jj^l- (b) Turbulence 
amplitude assumed, (c) Poloidal velocity profile (solid line), and velocity estimated 
using Doppler reflectometry signal frequency spectrum shift: • — using antenna 
focusing, o — without focusing. 

(p ~ 1 cm) to provide condition (fT2|l to be satisfied. Besides that we take into account 
the turbulence suppression in the barrier region (see figure Efb)). Despite the fact that 
density profile in the barrier region is favorable for the diagnostics, antenna focusing 
makes the spectrum shift more adequate to the behavior of plasma velocity in the cut-off. 

5. Discussion 

First of all we discuss the frameworks of approximations used. In this paper we 
consider Doppler refiectometry in slab plasma geometry. As it was mentioned in 
section ^ this model is reliable for large vertically elongated plasma. The effects of 
cylindricalness become important when the probing beam width in the cut-off vicinity 
is comparable with the cut-off surface curvature radius which takes place in small 
toroidal devices or in case of probing of plasma central regions. These cases of essentially 
cylindrical plasma geometry were considered numerically in ^3] and analytically (for 
specific plasma density profiles) in ^3]. In these papers rather obvious conclusion was 
obtained that plasma poloidal curvature enhances the diagnostics sensitivity to the 
fiuctuations with high poloidal wavenumbers. Another cylindrical geometry effect was 
investigated numerically in where strong influence of the plasma cylindricalness on 
the diagnostics locality was demonstrated. 

Essential can be the plasma poloidal curvature influence on the focused antenna 
beam. In the present paper it was neglected, which is correct, when the focus radial shift 
due to refraction associated with the plasma curvature is less than radial fluctuation 
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Here Aq is the ray trajectory length from the antenna to cut-off, which was defined above, 
Xc is the distance between the antenna and the cut-off and tq is the plasma radius in 
poloidal plane. Roughly this criterion can be approximated as x^/ro < icx- When it fails 
to hold to provide the focusing to the cut-off the probing wavefront curvature radius 
should be calculated taking into account plasma poloidal curvature. All the mentioned 
effects will be taken into account in separate paper of the authors. 

Additionally it should be noted that 2D plasma geometry effects like cylindricalness 
are not a priory significant in reflectometry (see e.g. [01 El EI, where it was 
demonstrated both analytically [HI and numerically |Tni E] that 2D theory predictions 
for radial correlation reflectometry are similar to conclusions of simple ID model). 

The second essential assumption made in the paper is associated with linear 
character of the scattering signal formation. According to this approximation is 
correct when following criterion is satisfied 



Here 6n is r.m.s amplitude of the turbulence and ric denotes the density in the cut- 
off. It can be seen that (jHUj) can be violated in case of high fluctuation amplitude or 
long trajectory length. In this case probing wave multi-scattering should be taken into 
account. This situation can be treated analytically by the procedure using in [TH] and it 
will be done in the paper by the authors which is now under preparation for submission. 

For diagnostics results interpretation it is important to distinguish linear and 
nonlinear situation. To do that experimentally one can use additional acquisition 
antenna (antenna 2 in flgure which receives the wave reflected from the cut- 
off. If the specular component persists in the frequency spectrum measured by this 
additional antenna the distortions of the probing wave in propagation are weak and 
we deal with linear situation of single-scattering. In the opposite case when the line 
at probing frequency is lost in the broadened reflection spectrum these distortions 
due to propagation in turbulent plasma lead to extinction of the specular component, 
which indicates transition to nonlinear small-angle multi-scattering regime. This way of 
experimental conflrmation seems to be reliable but needs additional access to plasma. If 
it is impossible some information can be obtained from the form of Doppler reflectometry 
spectrum received by antenna 1 in flgure ^ If the frequency spectrum width is consistent 
with estimation made based on the antenna pattern width (5uj ~ 2\/2v/p) one can 
conclude that the signal most likely resulted from single-scattering. 

The scattering efficiency introduced and analytically obtained dH)), in the 
paper reveals main similarities and differences of conventional reflectometry and Doppler 
technique. In the both methods scattering signal is proportional to reversed square of 
probing wavenumber [0]. That provides the diagnostics with spatial localization, which 
can be rather poor in case of unfavorable density proflle. 




(36) 
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The back-scattering signal component formation is similar for tilted probing, which 
is performed in Doppler technique and normal probing which is specific for conventional 
method. In particular, in both cases the probing wave focusing to the cut-off, according 
to (EH), underlines the cut-off region contribution and should improve the diagnostics 
localization. It is worth mentioning that utilization of non-slab probing wave fronts 
have been already discussed in j^l dl 120] , however in these papers the front curvature 
was chosen close to the cut off surface curvature in order to improve the Doppler 
reflectometry wave number resolution |Sl ED! , or to reduce the 2D interference effects 
in the signals reflected from different parts of corrugated cut off surface jTH]. Finally it 
should be noted that 2D-focusing onto the cut-off surface, which can be easily realized in 
experiment, can provide even better localization than ID-focusing only possible within 
2D model considered in the present paper. 

Essential are the peculiarities of FS efficiency (j22I)- Both for Doppler and 
conventional reflectometry the cut-off contribution to the FS component is received by 
the most favorable part of the antenna diagram ky = —q/2. However the suppression of 
the signal coming from the plasma volume in the case of conventional reflectometry takes 
place only for fluctuations satisfying condition gp ^ 1, at which the cut off contribution 
is suppressed as well. For longer poloidal scales > p the suppression is not efficient 
and additional localization is not possible. On contrary for Doppler reflectometry due 
to the tilted probing the FS component of the signal, coming from the plasma volume 
is suppressed providing the following condition is fulfilled 

It can be easily achieved by large enough antenna tilt angle or beam width, or by the 
focusing to the cut-off, thus making the FS contribution extremely localized to the 
cut-off. 

Finally we discuss the wavenumber resolution of the Doppler reflectometry. The 
scattering efficiency obtained (j^D), (PT|) demonstrates that diagnostics possess practically 
no radial wavenumber resolution due to the fact that scattering signal is an integral over 
radial wavenumbers and small x are pronounced with weight function 1/x. Poloidal 
wavenumber resolution can be easily estimated and is determined by antenna pattern 
width. For BS contribution it can be represented as 

oq ~ 

P 

FS poloidal wavenumber resolution can be \/2 times worse. 



6. Conclusion 



In the present paper the Doppler reflectometry spatial and wavenumber resolution is 
analyzed in the framework of the linear Born approximation in slab plasma model. The 
results obtained provide realistic description of Doppler reflectometry experiment in 
large elongated plasma at low level of density perturbation. 
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Explicit expression for the backscattering spectrum is obtained in terms 
of wavenumber and frequency spectra of turbulence assumed to be statistically 
inhomogeneous in radial direction. The treatment is performed for arbitrary density 
profile and antenna pattern taking into account diffraction effects. In agreement with jS] 
it is demonstrated that the signal consists of contributions of back and forward scattering 
in radial direction, which take place both before and after the reflection of the probing 
wave in the turning point. Similar to the traditional fluctuation reflectometry jH] 
the scattering efficiency for both back and forward scattering is shown to be inverse 
proportional to the square of radial wavenumber of the probing wave at the fluctuation 
location thus making the spatial resolution of diagnostics sensitive to the density profile. 
It is shown that additional localization is provided in case of forward scattering in the 
radial direction by the antenna diagram and in case of backscattering by the fact that 
the turbulence spectrum is suppressed at high radial wavenumbers. The improvement 
of the diagnostics locality by probing beam focusing onto the cut-off surface is proposed 
described as well. 

It is demonstrated that analytical expressions obtained can be easily used for 
fast interpretation of Doppler reflectometry data in particular for estimation of this 
diagnostics locality. They can as well serve for benchmarking and testing of full wave 
numerical codes developed for interpretation of conventional fluctuation reflectometry 
data [ini nil CSl Cni CZl EI] , should the authors consider the application of these codes 
to the field of Doppler reflectometry. 

The magnetic surfaces curvature, if important, can be accounted for within the same 
theoretical approach applied to cylindrical plasma geometry. The nonlinear effects in 
Doppler reflectometry spectra formation becoming significant in large devices and at 
high density perturbation level can be described within WKB approximation in the way 
similar to one used in |18j . 
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Appendix A. 

Here we estimate additional contribution to the forward scattering efficiency ()22|). which 
can be provided by singularities of the term |A(x, q)f / |x| in |T7j) . This term possesses 
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four branching points which in ray tracing case are determined by equations 



Xi_4 



"An -7^ 



cZ/* A* 



1/2 



1 ± 



2c(Ao - 7^) 



where L*, A^, are taken in the scattering point a;* (xi_4) corresponding to Xi_4. The 
branching point contribution is especially large in the cut-off vicinity where k'^{x) = 
oj"^ /c^{xc — x)/ and 

A(a;) = t^k{x) 



UJ 

The real part of right branching points = 9ftxi^4 in this case coincides with the 
stationary phase point x* = x^ _^ which gives the main contribution to the FS 
efficiency. The analysis in this case becomes complicated and inaccurate in WKB 
approximation. More rigorous approach to this case will be developed below in sectional 
taking into account that close to the cut-off, where density profile can be supposed linear, 
accurate solutions of equation are available. 

In general case these branching points are situated far from the stationary point, 
so that their contribution to the integral ()17j) is a quickly oscillating function of x and 
q and therefore is negligible. This is easy to show already in the case of profile slightly 
different from linear when branching point x^ is not so far from the stationary one to 
allow us to decompose 

x*(x6) ^ a + x'^(x{, - X*) 

To estimate the contribution of branching point we perform the integration over 
X = X — X;, within vicinity of the branching point taking into account that 



|A(x,g)r 



X 



c(Ao - 7^) 
The phase in fll7|l takes the form 

i[x-x^{x,q)f iuk'^{x)ll 



(2c(Ao - 7^) 



Ujp^Xij 



-1/2 



2x' 



ujK 



2cL^k{x) 



q' - 2ikXx 



c(Ao - 7^) 

This representation allows us to perform the integration over x and estimate the integral 
over q in ()17p. The ratio between the contributions of the branching and stationary 
points takes the form 



Sbjx) 



P 



UJ 



c(Ao - 7^) 



1 



ujA^ 



2cL^k{x) 



-1/2 



The factor on right-hand side of this equation can be estimated as 



c^A* 



2cL^k{x) 



-1/2 



UJ 



ck{x) 
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a 



where a = \Ll/ne{xc) ■d^ne(x)/dx^ 



characterizes the nonlinearity of the density 



profile. This leads to the condition (j24j) . when contribution in question can be neglected. 
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